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ABSTRACT 

We report results to date on the development of an autonomous, on board, 
near real time spacecraft attitude estimation technique. The approach will use 
CCD based star sensors to digitize relative star positions. Three microcomputers 
are envisioned, configured in parallel, to (i) determine star image centroids and 
delete spurious images, (ii) identify measured stars with stars in an on board 
catalog and determine discrete attitude estimates, (iii) integrate gyro rate meas- 
urements and determine optimal real time attitude estimates for use in the con- 
trol system and for feedback to the star identification algorithm. 

We present improved algorithms for the star identification. The discrete attitude 
estimation algorithm recovers thermally varying interlock angles between two 
star sensors. The optimed state estimation process recovers rate gyro biases in 
addition to real time attitude estimates. 
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Discussion of Slide Material 


Figure 1. 

For this technique, two star sensors are fixed to the spacecraft and 
to each other with interlock angles of, say, 60'’-90®. Both fields of 
view (FOV) are sampled at the same instant in time and the star data are 
used to determine the vehicle attitude. We note that each field can be 
used independently to yield attitude but if both fields of view are 
functioning, the data are combined to yield an orientation with higher 
precision and to recover interlock angles as well. We expect to attain a 
pointing accuracy of at least ^ 10 arc-sec. and probably S 5 arc-sec. 

Figure 2. 

This figure shows schematically the parallel processing nature of 
this technique. Process A will cycle independently of Processes B and 
C since it will likely be the most rapid. Process B will obtain image 
coordinates from Process A at the beginning of each cycle and at the 
same time obtain estimated orientation angles from the Runge-Kutta 
integrator section of Process C. The rion time of Process B will probably 
be the longest of the three but will also be quite variable. The output 
of Process B consists of three orientation angles and the covariance 
matrix obtained from the least-squares correction. These are used by 
Process C along with the estimated angles. We propose a Kalman filter 
algorithm to determine the optimal state estimate at the time associated 
with the current data from Process A. Process C will also update the 
estimated rate gyro biases. These are then used for the next Runge- 
Kutta integration. Process C will likely have a relatively fixed time 
budget . 
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Figure 3. 

Figure 3 displays a stammary of functions performed by Process A. 

CCD integration time is typically 0.1 - 1.0 sec. and readout time is < 0.1 
sec. The primary function of Process A is to determine image centroids. 
Star images are defocused slightly to spread images over 3 x 3 or 4 x 4 
pixels. The response of each pixel is used to interpolate positions to 
an accuracy of, say, one- tenth of a pixel. 

A second important function is to detect and delete spurious images . 
Images which are moving relative to other stars will be deleted. Since 
the cycle time should be quite short we can track images over several 
frames before Process B requests new data. 

The output of Process A will consist of image coordinates and 
brightness at the latest epoch, for all stars with intensity above some 
threshold. Each FOV is treated separately and passed to Process B. 

Figure 4. 

Process B is the most challenging process. Input data consist of 
image coordinates from each FOV and estimated attitude angles for the 
time associated with the image measurements. An on-board star catalog 
is accessed to yield coordinates for all catalog stars within some 
uncertainty border around the estimated FOV. These are sorted according 
to angle from the estimated boresight of the star sensor. 

Angles between measured stars are compared with angles between 
catalog stars until a tentative match occurs . The three orientation 
angles are adjusted to bring computed image coordinates of catalog stars 
into agreement (in a least- square s sense) with the image coordinates 
from the FOV. We next test for additional star matches and repeat the 
adjustment. 
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The above steps are performed for both FOVs. If both sensors are 
functioning, the data are then combined to refine the attitude angles 
and the three interlock angles between FOVl and F0V2. 

Output of Process B is three orientation angles at the time associ- 
ated with the image measurements. 

Figure 5. 

The image placement on a CCD star sensor is a function of six 
parameters: spacecraft orientation angles star coordinates (a, 

6), and the sensor focal length, f. The colinearity equations of Figure 
5 express this relationship for (x,y) placement. 

To simplify the computation we chose to store direction cosines for 
the star coordinates rather than (a, 6). Note that once [C] is computed 
for an estimated orientation, image coordinates are simple algebraic 
equations . 

Figure 6. 

In process B, the cosines of interstar angles are computed and 

compared. Computation involves only algebraic expressions containing 

•>> 

no trignometric functions. The star catalog contains direction cosines 
for each star rather than right ascension and declination (a, 5). Al- 
though more memory is required by the catalog, the speed gain by elimi- 
nating many sine and cosine calculations justifies the storage of 
direction cosines. 

Figure 7. 

Once a tentative star pair cosine from the catalog matches a 
measured pair cosine, the estimated orientation is iteratively adjusted 
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by least-squares differential correction. This scheine is used for both 
fields and the combined solution. In the latter case, three interlock 
angles are solved for as well as the three orientation angles. 

Figure 8. 

Figure 8 displays the stars accessed from the catalog for a typical 
case. The dashed line is the estimated FOV and the solid line represents 
the true FOV. The boxed stars are the stars that provided the first 
match with the measured stars. 

The sky has been divided into 256 non-overlapping cells, all nearly 
the same size. One or more cells are accessed for stars near the bore- 
sight. Only those stars that are within an uncertainty border (~11°) 
are used in the pairing logic. In this example, three cells were accessed 
giving 73 stars of which 38 were used in pairing logic. A match was 
found on the 38th attempt and, after correction, all four measured stars 
were matched with catalog stars. 

Only one cell, containing 24 stars, was needed for F0V2. A match 
was found on the fourth attempt and l’' additional star of the 5 measured 
stars was identified in the catalog after correction. The remaining 2 
measured stars were not contained in the catalog. 

Table 1. 

In this table we list the angle corrections obtained by a typical 
least-squares solution in Process B, The three attitude angles have 
been perturbed by approximately three degrees and the interlock angles 
by 1 arc-minute. 

As can be seen from the table, most of the attitude corrections are 
obtained by the FOVl solution. The attitude is eorrected again by the 
F0V2 solution due to poor interlock angles . Finally, the combined 
solution recovers both the attitude and interlock angles . 
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Figure 9. 


Figure 9 presents an outline of the main function of Process C. The 
actual algorithms to be used have not been finalized pending evaluation 
of Process B execution time in a microcomputer system. Short run times 
for Process B will allow less complexity in Process C algorithms. 

Presently, there are two primary tasks for Process C: 

1) Integrate the equations of motion ahead to the time of Process 
A output to B, This will use the last output of Process C for 
the initial attitude state and gyro rates will have biases 
subtracted. 

2) Use Extended Kalman Filter formulation to obtain the optimal 
best estimate of the spacecraft attitude at the instant of 
Process A output. Optimal estimates of gyro biases are also 
calculated. This algoritj|^m requires the estimated angles from 
1) above, the discrete attitude arid associated covariance matrix 
obtained from the Process B solution, and the state covariance 
matrix from the last Kalman filter update. The best estimates 
of angles and biases are then used for subsequent integration 

of the motion equations. 

Figure 10. 

This figure demonstrates the behavior of the Kalman filter on simulated 
data. The plotted data show the deviations of the first attitude angle 
((j)) from the true orientation. Attitude from Process B is simulated by 
perturbing the true orientation angles (standard deviation of 8.8 arc-sec). 

The gyro rates are perturbed from a constant value (-4 x 10 rad. /sec.; 

—6 

Standard deviation of 1.5 x 10 rad. /sec.) each time they are needed and 

—6 

biases are set to 4.8 x 10 rad. /sec. The initial estimated state is off 
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10° in each angle from the truth (very conservative) and biases are set 
to zero. 

The attitude state is integrated ahead in 1 minute intervals and 
the deviation from the truth calculated (shown as Simulated data 

from Process B then becomes available and the Kalman filter applied. 

The deviations of the filtered results from the truth is denoted by 
crosses (*'+”) . Due to noisy gyro rate data and the 1 minute integration 
time, the filtered attitude nearly equals the Process B attitude. 

The dashes (”~”) are the estimated deviations as calculated from 
the updated covariance matrix after each pass through the Kalman filter. 

Note that after only three Kalman filter applications the attitude 
and biases have been recovered, and the steady state attitude errors are 
less than 10 arc-seconds in all cases and usually less than 5 arc- 

s 

seconds. The, one minute interval between star-determined orientation 
estimates (from Process B) will in practice be variable and typically 
less than one minute; the precision of these simulated converged esti- 
mates is therefore felt to be conservative. 
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FIGDKE 1 DVASTAR An electro-optical/software system capable 
of real ti me readout of digitized star coordinates, 
and ultimately, autonomous, near-real time star 
pattern recognition and attitude determination. 
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Figure 3: Major tasks of Process A. 
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Figure 4: Major tasks of Process 
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Figure 7 : Method of least squares for estimated attitude 

angle correction. are orientation 

angles and (x,y) are image coordinates. 
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Figure 8: Boundaries of the true field of view (solid rectangle) and 

estimated field of view (dashed rectangle) . Center of 
each FOV is denoted by The three large sectors are 

the three cells accessed by Process B. The pair of boxed 
stars is the first pair to match the measured stars. FOV 
is approximately 7* x 9*. 
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Figure 1.0; Behavior of Kalman filter for first 




Figure 11: Concluding Remarks 


A new spacecraft attitude estimation approach is under development. 
Key features are: 

*Electro-optical star sensors with no moving parts 
*0n board, near real time attitude estimation 
*Sub ten arc-second precision appears feasible. 


A three parallel process division of labor is proposed: 


A - Star-sensing/image processing 

B - Star identif ication/Discrete attitude estimation 
C - Motion integration/Attitude estimation 


The software for Processes B and C is currently being implemented on 
microcomputers. The algorithms will be tested, evaluated and modified 
to maximize data throughput and accuracy. 

Area of major effort and/or problems: 

*Limited precision and speed obtained with on-board microcomputers 
*Data acquisition rate vs. Process B run time 
*Process B run time vs. Process C complexity. 
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